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We introduce a model to describe columnar phases of chiral discotic liquid phases in which the 
normals to disc-like molecules are constrained to lie parallel to columnar axes. The model includes 
separate chiral interactions favoring, respectively, relative twist of chiral molecules along the axes of 
the columns and twist of the two-dimensional columnar lattice. It also includes a coupling between 
the lattice and the orientation of the discotic molecules. We discuss the instability of the aligned 
hexagonal lattice phase to the formation of a soliton lattice in which molecules twist within their 
columns without affecting the lattice and to the formation of a moire phase consisting of a periodic 
array of twist grain boundaries perpendicular to the columns. 



PACS numbers 61.30.Cz,61.30.Jf,61.72.Bb 
I. INTRODUCTION 

Chirality gives rise to a rich variety of modulated 
equilibrium phases in liquid crystals B, including the 
cholesteric, smectic-C*, TGB, and blue phases. To date 
most theoretical and experimental work has focussed on 
systems composed of molecules that adopt rod-like con- 
formations that favor the formation of uniaxial smectic 
phases. Discotic chiral molecules favoring the formation 
of columnar phase have, however, been synthesized [^| -^). 
They produce discotic cholesteric phases [|J and colum- 
nar phases with chiral structure. They also exhibit inter- 
esting ferroelectric properties 

Chirality favors twisted structures. In chiral smectic 
phases, twist can be expelled altogether in the smectic- A 
phase, it can appear as molecular twist in the smectic-C* 
phase, or it can appear as layer twist in the TGB phases 
||. Columnar phases can exhibit the analog of all of 
these phases and some phases with no analog in smec- 
tic systems. Some possible columnar discotic phases are 
shown in Fig. [I]. If chiral forces are sufficiently weak, the 
lattice structure of the columnar phase can simply expel 
twist. If the coupling between molecular chirality and the 
columnar lattice is sufficiently strong Q , chirality can in 
principle induce a tilt-grain-boundary phase, analogous 
to the TGB phase in smectics, with rotation axis perpen- 
dicular to the columns, or a moire phase with rotation 
axis parallel to the columns. In the former phase, there 
is a periodic lattice of grain boundaries separating ro- 
tated regions of aligned columns. In the latter, there is 
a periodic lattice of grain boundaries perpendicular to 
the columns across which the hexagonal columnar lat- 
tice undergoes discrete rotations. In columnar phases, 
molecules can twist within their columns without de- 
stroying the lattice structure. Such behavior is possi- 
ble in columnar systems because each molecule sits in a 
fairly symmetric environment, and it can rotate without 
greatly disturbing its neighbors. In smectic phases, the 
rotation of a molecule within a layer would lead to enor- 
mous disruptions and would be energetically very costly. 
Phases with spontaneous chiral symmetry breaking in 
which propeller-like molecules rotate in different direc- 



tions in different columns have been observed experimen- 
tally (8) and modeled theoretically M. In chiral systems, 
the direction of preferred rotation in all columns is set by 
the underlying molecular chirality. The pattern of molec- 
ular rotation in these phases can be described as a soliton 
lattice, and we will refer to them as soliton phases. 





(b) 





(c) (d) 

FIG. 1. Some chiral discotic columnar phases: (a) A uni- 
form orientationally ordered phase in which there is not twist. 
This twist has been expelled as in the smectic- A phase of chi- 
ral molecules, (b) An orientationally disordered or "plastic" 
columnar phase. The is no long-range orientational order of 
the molecules, but the hexagonal columnar structure remains 
intact, (c) A cholesteric-like phase in which molecules ro- 
tate in a helical fashion within each column while maintain- 
ing phase coherence between columns, (d) The moire phase 
in which the columns themselves rotate about the 2 axis. Ro- 
tations occur in discrete jumps about twist grain boundaries. 
(From Ref. §) 
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In this paper, we will introduce a fairly general model 
for chiral columnar phase that has soliton and moire 
phases in addition to phases with expelled twist. To 
keep the model as simple as possible, we restrict the 
normal to the disc-like molecules to lie parallel to the 
columnar axes. We will, therefore, not be able to discuss 
the tilt-grain-boundary phase or ferroelectric chiral dis- 
cotics. In a future publication, we will study a variety of 
phases that can result with this constraint is relaxed. Our 
model begins with interacting chiral molecules in a dis- 
cotic columnar phase. Its Hamiltonian has an elastic part 
associated with distortions of the lattice. It has terms 
favoring parallel alignment of neighboring molecules and 
chiral terms favoring both molecular and lattice rotation. 
It also has crystal-field terms coupling molecular orien- 
tation to the columnar lattice and defining a twist pene- 
tration depth Xg analogous to that of smectic systems. If 
lattice distortions are prohibited, our model is essentially 
identical to that studied by the Sherbrooke group || . We 
assume that the low-temperature, weak chirality, ground 
state of our Hamiltonian is the ordered "ferromagneti- 
cally" aligned state shown in Fig. [l]a, though other large 
unit-cell states are possible || . As chirality is increased, 
the aligned state can become unstable either to a soliton 
lattice if the crystal-field coupling is weak or to a moire 
phase if it is strong. Our primary aim will be to deter- 
mine the critical chirality for these two instabilities. Our 
model, however, has the potential for more complicated 
phases and very complex phase diagrams. For example, 
a mixed moire-soliton phase in which there is molecu- 
lar twist relative to the lattice between grain boundaries 
could exist. Or the soliton phase could melt altogether 
to a "plastic" discotic phase with no orientational long- 
range order (Fig. 0b). This phase could them become 
unstable with respect to the formation of a moire phase. 

Though our model is motivated by chiral discotic liquid 
crystals, it can, with proper interpretation, be applied to 
aligned chiral polymers such as DNA JTo| . Aligned poly- 
mers can form a hexagonal lattice perpendicular to the 
direction of alignment. A polymer like DNA is a tightly 
wound double helix. This structure makes it unlikely for 
it to form an oricntationally ordered phase analogous to 
the discotic ground state shown in Fig. |l|a. Rather the 
phases of the molecular helices on neighboring molecules 
will be random: there will be no long-range order in 
molecular orientations perpendicular to the direction of 
polymeric alignment. The ground state will be equiva- 
lent to the disordered states shown in Fig. [l]b, though 
each molecule will have an average twist. Thus, there is 
no analog of the soliton phase in DNA, and the moire 
phase is driven by chiral terms favoring lattice rotation 
rather than molecular rotation. 

The outline of this paper is as follows. In Sec. II, 
we define the model and discuss its continuum limt. In 
Sec. Ill we discuss instabilities toward the moire phase. 
This analysis differs from that of reference Q] because 
of the finite twist penetration depth. Our results are, 
however, almost identical to those of that reference. In 



Sec. IV, we discuss the soliton phase and arrive at a 
criterion determining whether the soliton or the moire 
phase will form. 

II. THE MODEL 

A. The Lattice Model 

Columnar forming chiral molecules can come in many 
forms. We will limit our discussion to molecules such as 
that shown in Fig. || with C3 symmetry. This molecule 
is similar to some that have recently been synthesized Q 
and to those spontaneously formed in the experiments 
of Heiney Q . Since we are interested principally in the 
competition between soliton and moire phases, we will as- 
sume that the molecular normal always aligns along the 
columnar direction. We will, therefore, not be able to dis- 
cuss the formation of the tilt-grain-boundary phase. We 
use a mixed continuum-lattice description of the colum- 
nar phase. The discotic columns form a two-dimensional 
hexagonal lattice with lattice parameter a. They are la- 
beled by an index 1. Distance parallel to the columns is 
specified by the continuous coordinate z. The orienta- 
tion of molecules at position z in column 1 is specified 
by the angle Q\(z). The column coordinates are given by 
Xi(z) = Ri + u(l, z) , where Ri is the equilibrium two- 
dimensional lattice coordinate and u(l, z) is displacement 
from equilibrium that can depend on both z and 1. 




FIG. 2. Schematic representation of a chiral discotic 
molecules with C3 symmetry. The three vanes coming off 
the hexagonal core are all tilted in the same direction relative 
to the normal like the blades of an airplane propeller. 

The Hamiltonian for our lattice model can be divided 
into an elastic part TL C \, an angle part Tig, and a chiral 
part Ji . The elastic Hamiltonian is the standard one for 
a columnar structure. In the harmonic limit, it is 

Hci = - I dz Kjj (1, \')uj(\, z)uj(Y, z) 
1 J 1,1' 

+~ f dz]>>[c\ 2 u(l,z)] 2 . (1) 
J 1 

The first term in this expression, with Xy(l,l') an elas- 
tic constant tensor, is the familiar elastic energy of a 
two-dimensional harmonic lattice of columns. The sec- 
ond term, with k a bending rigidity, measures the energy 
of bending the columns. We assume that neighboring 



2 



molecules want to be parallel in the absence of chirality. 
Furthermore, there are couplings between lattice distor- 
tions and molecular rotation and a preferred orientation 
of the molecules relative to the lattice. All of these effects 
are incorporated into the model Hamiltonian 

He = -~A j dz^ cos[30i(z) - 39 v (z)} 
2 J i,i' 

-]-B j dz^cos[66»i(z) -60i. r (z)] 
2 J i,i' 

~C j dz^Gw[Wi{z) + W v {z)-&<k,v{z)], 
^ i,i' 

ld*\0M*)? (2) 
i J 

where the sum 1, 1' is over nearest neighbor columns in 
the lattice and 

,M*) = tan- ( f+<h,M (3) 
\fli + u x {\, z) - u x (V , z) J 

is the angle the bond between column 1 and 1' makes with 
the s-axis at z and a 1 ' 1 is the equilibrium lattice vector 
connecting those columns. The Hamiltonian He is invari- 
ant under 9\ — * #i+(27rn/3) for any interger n, as required 
by the C3 molecular symmetry. It is also invariant under 
rotations of the lattice by 2ir/6 and under simultaneous 
rotations of the molecules and the lattice through arbi- 
trary angles. When lattice distortions arc prohibited, He 
is very similar to that studied by the Sherbrooke group 

i-. 

Finally chiral interactions along a given column fa- 
vor molecular rotation, and chiral interactions between 
molecules in different columns favor lattice rotation. We 
introduce two chiral terms in our model to describe these 
effects: 



(4) 



K*9 = ~le j dz^2d z 6i(z) 
J 1 

ft-l = -!<$> / dzd z <k,Y- 
1 v J 



In the ordered phase of discotic systems, the dominant 
twist comes from Wg, and we may assume 75 3> 70- In 
polymeric systems or in the orientationally disordered 
phase, 76/ is zero, and any rotation is induced by Hj, 



B. The Continuum Limit 

When spatial variations are slow on the scale of the 
lattice spacing, we may expand the lattice Hamiltonian 
in gradients of lattice displacements u and angles 9 and 
replace the lattice sum by a continuum integral. To this 
end, we replace u(l, z) by u(x) and 6\{z) by 0(x), where 



x = (xi,2) with xj_ the coordinate perpendicular to z, 
and we set J dz J^i ~* a ' 1 J d 3 x. In this limit, the lattice 
angle 4>\ t y (z) can be expanded about its equilibrium value 
of 4>\ V {z) as 

= a~ 2 [a l x a\ l d l u y (x) - aj; 1 a'' 1 9,ii x (x)]. (5) 

The Hamiltonian He does not depend on the equilibrium 
angle 4>\\> because the latter is an integral multiple of 
27r/6. The average over nearest neighbor lattice sites of 
$<h,v{ z ) IS the hexatic angle <^>6( x ) : 



4>6 (x) 



2 ^13 v l U,j 



(6) 



r 



where tij is the anti-symmetric two-dimensional tensor 
with i and j running over x and y. In addition, 



53<ty&'(*)] = 6$ + 3u 



tjUij 4 U U 



(7) 



where itjj — (diUj + djUi)/2 is the symmetrized strain. 
With this information, we can express the continuum- 
limit Hamiltonian as a sum of terms: H = H c \ + He + 
H g + H* . The elastic Hamiltonian 

We! = \ J d 3 x[Xul + 2fiu ljUlJ + K(d 2 u)% (8) 

is the standard continuum elastic Hamiltonian for a 
hexagonal columnar system. Here K — k/vq, where 
v = \fia 2 /A is the area of an hexagonal unit cell, and the 
Lame coefficients A and fi are determined by the contin- 
uum limit of Kij(l, 1') and by the the parameters B and 
C in He through the second order expansion in 5(()\ } \>. 
The angle Hamiltonian is simply that of an anisotropic 
xy model: 



He 



(9) 



where K e l = n e /v , = 9 (A + C)/2, and V± = 

(d x , d y ,0). The coupling term H g is in the small 6 — (pe 
limit is 



H g = \g I d A x 



3 ) 2 



(10) 



where g = 108(S+C7) / a 2 and where the final form is valid 
in the small 9 limit. Finally, the chiral energy becomes 



H* 



-7e / d 3 xd z 9-j,p / d 3 xd z (p i: . 



(11) 



where je = Je/ v o and 74, = j^/vq. 

The continuum Hamiltonian has a structure imposed 
by rotational invariance: the coupling term \g{9 — 
^eijdiUj) 2 is invariant under simultaneous rotations of 
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the lattice and molecular orientations. It is the ana- 
log for columnar systems of the invariant coupling Jll]] 
(<5n + Vw) 2 of smectic-A liquid crystals, where <5n is the 
deviation of the Frank director form its equilibrium ori- 
entation. The columnar phase, like the smectic-A phase 
tends to expel molecular twist. In the ground state, 
molecules align along preferred crystal axes with 9 = 0. 
A twist in 9 relative to the lattice at z = will decay to 
zero in twist penetration depths 



A 



II 



Kjj/g, Xg = JKfr/g 



which tend to zero in this strong coupling, g 
Another important length is 



(12) 



oo limit. 



(13) 



giving the length scale over which bend deformations 
heal. 

Two limits of the continuum Hamiltonian deserve spe- 
cial attention. The first is that in which there is no long- 
range orientational order. In this limit, 9 looses its mean- 
ing, and only terms involving u remain. The Hamilto- 
nian becomes Ti e \ + H.^ with, in particular, g = 0. This 
is the limit studied in fij. The second is the strong- 
coupling g — > oo limit in which 9 and 0g are forced to 
be equal. This leads to H* = — (70 + 70) J d 3 xd z 4>e and 

He = ±J d 3 x[Kj-{V xeijdwj) 2 + K^d^diUj) 2 }. Thus, 
except for a (d 2 u) 2 correction, the form of the Hamil- 
tonian with g = 00 is identical to that for g = 0. Wc 
should expect, therefore, that the critical values of chi- 
ral couplings leading to the moire phase will have simi- 
lar forms but differ in magnitude in the two limits and 
that there will be a smooth, non-singular interpolation 
between these limits as a function of g. 



III. THE SOLITON PHASE 

If the coupling g between molecular and lattice ro- 
tations is weak, molecules will be able to rotate rela- 
tive to the fixed lattice. We can describe this situation 
by a Hamiltonian depending only on 9 and not on u: 



H = H e +H q +'H*n, where 



Tig = — 5^5 / d d xcos69. 



(14) 



This Hamiltonian is a chiral Sine-Gordon model that 
is equivalent to the continuum limit of the Frenkel- 
Kontorowa model |L2j and that describes a cholesteric 
liquid crystal in an external field. The chiral term Tig 
favors twist that the angle elastic term Hg oposes. Twist 
sets in when the coupling g exceeds the critical value 
necessary to nucleate a single soliton, which produces a 
rotation of 9 through an angle of ir/3 from one end of the 
sample to another. The energy per unit area of a single 
soliton is 



(15) 



The chiral energy gain arising from a single soliton is 
Hg = —A(tt/3)^b, where A is the area. Thus the total 
energy per unit area of a single soliton is 

/. = W - (*79/3)]. (16) 
The critical chiral coupling constant is, therefore, 
. 2 



(17) 



For 7e > 7^, there will be a soliton lattice in 9 with a 
lattice spacing that decreases with increasing 79. This 
is a helical state with a pitch equal to the soliton lattice 
spacing. 



IV. THE MOIRE PHASE 

The moire phase consists of a periodic array of twist 
grain boundaries perpendicular to the columnar axis 
across which the orientation of the hexagonal lattice ro- 
tates in discrete jumps. The twist grain boundary is 
a honeycomb lattice of screw dislocations. This phase 
forms when the energy cost of creating a low-angle grain 
boundary is just counterbalanced by the twist energy 
gain arising from TL* . Let lb be the distance between 
grain boundaries and l d be the length of the side of the 
hexagonal unit cell of the honeycomb dislocation lat- 
tice (see Fig. ||). The total length of screw disloca- 
tions in a lattice of N cells is L = 3Nld- The total 
area of this lattice is A = 3Nl 2 /2\/3. In the limit of 
large lattice spacing (i.e., large l b and Id), interactions 
between dislocations can be neglected, and the energy 
per unit volume of an array of low-angle grain bound- 
aries is f GB = e d L/(Al b ) = 2e d /(l d y/3l b ), where e d is the 
energy per unit length of a dislocation. Far from the 
boundary, 9 = (()§, and both 9 and ipe undergo the same 
jump A9 across the boundary. This jump was calculated 
in Ref. 0. It is b/l d V3 where b = d is the magnitude 
of the Burgers vector. The chiral energy of an array of 
grain N boundaries is thus 



NA 



V- 



jb 

l d V3 V^lbld 



(18) 



where 7 = jg + 70 and V = NAlb is the volume. Thus, 
the energy per unit volume of the moire phase when dis- 
locations interactions are ignored is 



fni 



1 



■(2e d - 7 6). 



(19) 



This energy becomes negative, and the moire state be- 
comes energetically preferable to the ordered phase when 
7 > 7™, where 



4 



7" 



2e d /b. 



(20) 



This result does not depend on the particular type of 
dislocation lattice formed in the grain boundary. For 
example, if the grain boundary consists of identical or- 
thogonal grids of dislocations with separation l d between 
dislocations then L — 2l d , A — 1%, and A9 — b/l d so that 
f m = (2e d - -fb)/(l d l b ), again producing 7" 1 = 2e d /b. 



(25) 



(26) 



where tpi = i(qi<J — e^qjit) / 'q\ defines the longitudinal 
part, which is obtained by solving Eq. (24). Substituting 
Eq. ( p5| ) into Eq. ( pjj ) and using q^i = ia and tuqiipi = 
wr, we obtain 



FIG. 3. Honeycomb lattice of screw dislocations that com- 
prise twist grain boundaries in the moire phase 



1 



(2/i + \)q\ + Kq\ 
2^ - Kq\ 



-£imqiqjq z a ri 



1 



(/3(g) + fi)ql + Kqi 

x [-P(<l)q±(q±qna Z n{q) + q±q z a nn {q) 



(27) 



(28) 



H + Kq 2 z fi . 
2 qzqiqieije im a mj ^qiqiq P ^ijej pu Oi U i[q) 



where 



The ordered phase becomes unstable with respect to 
the formation of the moire phase when 7 exceeds 7 m . Our 
task, therefore, is to calculate the energy per unit length 
of a screw dislocation e d . We follow closely the procedure 
of Ref. jl3| appropriately generalized to include 9 as an 
independent variable. We introduce u> 7 i, which is equal 
to d-yUi away from defects, where i = 1, 2 and 7 = x, y, z. 
We also introduce the dislocation density, 



a 7 i(x) 



E 



dst nn (s)b n J i3) (R n (s) -x), (21) 



where R n (s) is the position vector of dislocation n with 
Burgers vector b„ as a function of its arclength s and 
tn(s) = dR n (s) / ds is its unit tangent vector. The con- 
dition <f> gr dui — bi that the integral of the changes in u 
around a contour dT enclosing a dislocation with Burgers 
vector b be equal to b then implies the constraint, 



^fx^Oi/W^i — Q^(x), 



(22) 



on Wyi. To find the elastic energy associated with dislo- 
cations, we need to minimize the energy Ji = 7i c \ + Tig 
subject to this constraint. Minimizing Ti with respect to 
variations in 9 and u, we obtain 



SH 



^ = -KfVi_e - ffiivjj* + 9(6 - y ijWij ) = 0, (23) 



— = -n{d iWij 



XdjWu 



K3lw ZJ 



+ -(eijdit 



-eijCij/diiWi/jt) = 0. 



After Fourier transforming, we can solve Eqs. 
( p2| ) for 9 and w 7i f. 



(24) 
and 



1 g(K]}qf\+Kj-<? ± ) 
±g + K$q?,+Kj-qi 



(29) 



For a single screw dislocation aligned along the x-axis. 
a vi (q) = 2<Kb8 vx 6 ix 8{q x ). Using Eqs. © ©, and ©, 
we obtain 



Wy X (?) 

w zx {q) 



where 



-iq z axx{q) Kq 2 q 2 z 

q 2 TM)q 2 v + Kqi 
iq y axx(q) M(q)g 2 
q 2 H<i)q 2 y + Kqf 



Mq) = /•* + /%)■ 



(30) 



(31) 



Then, using Eq. ( p0[ ) in 7i, we obtain the energy per unit 
length of a dislocation, 



= -Kb 2 ( dq y dqz gMg) 

d 2 7 (2tt)2 7J(q) g 2 + ^^ 



(32) 



where the integrals over (/j, and q z have respective upper 
cutoffs of the inverse correlation lengths, ^J 1 and £,T . 
This expression is identical to that obtained in Ref. 
with 7i(q) replacing /i. 

To evaluate ed, it is convenient to express the integral 
in a unitless form: 



' ° lo dy l dZ y 2 d + f) + S^- 



z 2 {l + f) 



(33) 



where 



5 



eo = 



and 



1 + ot^z 2 + a\y 2 ' 



(34) 



(35) 



In the expressions for e d and /, we have introduced unit- 
less ratios, all of which are Ginzburg parameters mea- 
suring the ratio of a penetration depth to a coherence 
length. The parameter 



<5 = 



was introduced in Ref. j7j. The parameters 
Jl \± 



qi 



(36) 



(37) 



and 




011 



1 IKjt A 



(38) 



are Ginzburg parameters for twist penetration. Note that 
j3± and /3|| are independent of 3. 

Analytic evaluation of e c i is difficult except for 5 = 0. 
We will content ourselves with the <5 — > and 8 — > 00 
limits for g small and for g — 00. Since ded/dg is positive, 
6d will increase monotonically from the its value at g = 
to its value at g = 00. As g — ► 0, we find 



1. gr -> 0, 5 = 



e<2 



1 



1 K 1 ' 2 ^ 2 



16tt 



qi 



1 + ^ 

8^ 



(39) 
(40) 



and 
2. g -> 0, <5 



v/2 



e <T 



1 



£5 
16 



1 j^W / _3_£ 
8\/27r V?I V 16 At 



(41) 



For g — 00, we find 



1. .g = 00, <5 — > 

r-2 



1 



2. g = 00, <5 = 00 



(42) 



e d = ^o8-"\ J X dyy-WQ + P\ff" (43) 



These results reduce to those of reference (when a 
missing factor of 1/4 is added), ed increases smoothly 
and monotonically with g. Its value at g — 00 is finite 
and depends on the Ginzburg parameters /3j_ and (3\\. 
When /3m and /3± are both zero, / = 0, and e d (g) = 
€d(g = 0). When these quantities are much greater than 
unit, e d {g — 00) > e d {g = 0): 



e d (g = 00, S = 0) 
e d (g = 0,6=0) 
e d {g = 00, 8 = 00) 
e d (g = 0, 8 = 00) 



P\\ » 1 
P ± » 1. 



(44) 



Thus, large g and large angle elastic constants K g and 
Kjj- lead to large values of and suppress the formation 
of the moire phase. 



V. DISCUSSION AND REVIEW 

In this paper, we have developed a model for chiral 
discotic columnar liquid crystals, and we have investi- 
gated its instability toward the formation of two types 
of structurally chiral phases: the soliton phase and the 
moire phase. Chirality in our model, which restricts the 
average molecular normal to be along the columnar axis, 
gives rise to two kinds of chiral interactions, one tending 
to rotate molecules within a column and the other tend- 
ing to rotate the columns themselves. These two interac- 
tions are characterized by respective coupling strengths 
70 and 70. There is an energy cost associated with ro- 
tation of molecules relative to the lattice characterized 
by a coupling constant g. When g is small, rotation of 
molecules within columns with a fixed lattice structure 
is possible. This is the soliton phase that develops for 
molecular chiral coupling constant 70 greater than a crit- 
ical value 7| ~ yfg. When g is larger formation of the 
moire phase is favored for 70 + 7^ > 7™ where 7™ is a 
smooth function of g that is finite in the g — > 00 limit. 

We have focussed on the instabilities toward two pos- 
sible structurally chiral phases. A full phase diagram for 
the model and indeed for real chiral discotic systems can 
be quite complex with mixed soliton-moire phases. Addi- 
tional phases can occur if the constraint that the molecu- 
lar normal (the Frank director) be parallel to the colum- 
nar axis be relaxed. In particular, tilt-grain-boundary 
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phases and smectic-C*-like phase in which the director 
tilts relative to the columnar axis and rotates in a helical 
fashion along that axis can occur. These more complex 
phase are currently being investigated. 

We are extremely grateful for many conversations with 
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ger. This work primarily by the MRSEC Program of 
the National Science Foundation under Award Number 
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